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ABSTOACP 

The radiation balance of the earth-atmosE^re systen 
simulated by using the general circulation model (GCM) 
of the Laboratory for Atmospheric Sciences (GLAS) is 
examined in regards to its geographical distribution, 
zonally-averaged distribution, and global mean. Most 
of the main features of the radiation balance at the 
top of the atmosphere are reasonably simulated, with 
some differences in the detailed structure of the pat- 
terns and intensities for both simmer and winter in 
cornparison with values as derived from Nimbus and NCAA 
(National Oceanic and Ataoj^heric Administration) 
satellite observations. Both the capability and de- 
fects of the model are discussed. 

INTRODUCTION 

The GEAS GCM simulates the general features of climatology 
reasonably well (Halen et al . , 1978). aimatological elanents, 
such as the mass distribution, geopotential height, and oceanic 
rainfall rates generally are in reasonable agreement with ob- 
servations. Here, we examine the primary driving force, radia- 
ticxi balance, of the circulation system. The GIAS GCM has the 
feature that the important thermodynamic variables, such as cloud 
types and cloud heights, atmospheric and ground teiperatures, 
and specific hianidity, are interactive with model dynamics, 
which allows this study to be meaningful. Through a corplicated 
interactive process, in order for the simulated radiation 
balan<^ at the top of the atmosphere to be right, the kinematics, 
dynamics, and thermodynamics of the model have to be right. The 
geographical distribution of the earth-atmosphere system radia- 
tion balance from satellite measurements provides a valuable 
check on our diagnostic studies of the circulation system of 
the GLAS GCM. 

The winter (sumner) simulations, starting from real data 
for 1 January 1975 (15 May 1974), were integrated for 60 days 
( 105 days) , and the last 30-day mean was used for comparison 
with geographical distribution of the radiation balance as 
derived from Nimbus 3 measurenents, for the 14-day (16-day) 
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iT>ean» starting from 21 January 1970 (16 July 1969). The 
zonal mean of the radiation balance was cctnpared with both 
Nimbus 3 and NOAA SR measureiients (January 1975). 

It should be noted that the periods used for cctnparison 
are different since Nimbus 3 data was not available for the 
exact period of years for vrfiich the model had equilibrated. 

It should also be noted that the satellite measurement instru- 
ment is not calibrated on board fpr shortwave measurement; 
and it is a sun-synchronous orbit. 

The Geographical Distribution of the Radiation Balance of the 
Earth-Atmosphere System 

Global patterns of the radiation balance and its asso- 
ciated cai^jonents, absorbed solar radiation, and outgoing long- 
wave radiation, are shown for winter in Figs. 1, 2, and 3, 
respectively. Similar distributions for simmer are shown in 
Figs. 4, 5, and 6. For cornparison, valres as derived from 
Nimbus 3 measurement are displayed in Figs. la-6a, and values 
as simulated from the model are displayed in Figs. lb-6b. 

Overall, there is reasonably good agreement between the 
simulated and the observed radiation balance patterns, particu- 
larly for the major energy gain areas over the subtropical mini- 
mim cloudiness regions in both simmer hemisp^iere oceans, for 
local minima over the continents in sunmer hemisE^ieres, and for 
the zonal pattern of deficit in the winter hemispheres. For 
all of these features there is general agre«nent in patterns, 
which indicates that the GLAS GCM is capable of simulating the 
general climatology to the accuracy of ^ zeroth order. *1116 
differences are in the intensities and in the detailed struc- 
ture of the patterns. It is from these differences that we 
are able to diagnose some of the model defects. 

In winter, the simulated zero-balance isoline is around 
30N without a southward dip off the west coast of North America, 
probably because the model does not generate low stratus clouds, 
and with a relatively small southward shift of the zero-balance 
isoline over North Africa, probably because the model generates 
slightly more cloudiness over Africa especially on the equatoi^- 
ward side of North Africa. The simulated period is February, 
whereas the observed period is the last two weeks of January. 

Due to the sun's declination, there should be about 5* northward 
shift of the zero-balance isoline. The cause of the remaining 
latitudinal difference may be due to the following reasons. 
Generally speaking, the model overestimates total absorbed solar 
radiation in the winter hemisjAiere in ccmparison with observation 
(Fig. 2). In addition, the model systematically underestimates 
outgoing longwave radiation (Fig. 3). ^r example, around 20N, 
the simulated value is about .36 cal/cm /min, whereas the value 
as derived from observation is about .42 cal/cm2/)nin. 
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TOTAL ABSORBED SOLAR RADIATION (cal/cmVmin ) 


a. Observed (Jan.) 




Fig. 2. Total absorbed solar radiation in the earth-atmosphere 
system 

a, as derived from Nimbus 3 measurements during the 
period 21 January to 3 February 1970. 

b. as simualted from the GLAS GCM for February. 
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OUTGOING LONG WAVE RADIATION (cal/cmVmin ) 


a. Observed (Jan.) 



Fig. 3. Outgoing longwave radiation emitted from the earth-atmo 
sphere system to space 

a. as derived from Nimbus 3 measurements during the 
period 21 January to 3 February 1970. 

b. as simulated from the GLAS GCM for February. 
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RADIATION BALANCE (cal/cmVmin ) 




Total absorbed solar radiation absorbed ir the earth 
atmosphere system 

a. as derived from Nimbus 3 measurements during the 
period 16 to 21 July 1969. 

b. as simulated from the GLAS GCM for August. 



OUTGOING LONG WAVE RADIATION (cal/cmVmln ) 


a. Observed (July) 



b. Simulated (Aug.) 



Fig. 6. Outgoing longwave radiation emitted from the earth-atmo-* 
sphere system to space 

a. as derived from Nimbus 3 measurements during the 
period 16 to 21 July 1969. 

b. as simulated from the GLAS GCM for August. 
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The major energy gain area is over the Southern Hemisriiere 
with i^ima over the oceans in connection with subtropical 
cloudiness minima, and local minima over continents associated 
with a heavy ITCZ cloudiness over South America and South 
Africa around latitude 0 to 20S. The model reprtxiuces maxima 
over the oceans, because the model simulates subtropical anti- 
cyclcaies reasonably well, as well as local minima over the 
continents. 

As would be expected from seasonal change in solar heat- 
ing, the major energy gain areas are moved to the northern 
sunmer for July and August with maximim ever the oceans. One 
striking feature appearing in the sunmer c4>servations is a 
deficit in the energy budget over the Sahara and Saudi Arabia, 
•nie model produces a more localized minimun over the Sahara, 
in agreement with the wide-angle radiative balance for August 
measured with the Earth Radiation Budget (ERB) instriment in 
Nimbus 6, v^ich also does not show a negative deficit as in the 
earlier Nimbus 2 and 3 observations. Hiis suggests that part 
of the discrepancy may be due to the coarseness of resolution 
of the model. When we examine the separate cemponents of the 
radiation balance, we find that the deficit is due to high 
surface albedo and high outgoing longwave radiation. But the 
particular reason for the deficit to be over the northernmost 
part of Africa, betv^n 30N and 40N, can be found by examining 
Figs. 5a and 6a. Fig. 6a is a display of the outgoing longwave 
radiation. Ihe cloudiness associated with the ITCZ and stumer 
monsoon system are reasonably well revealed in this figure. 

Fig. 5a shows that the absorbed solar radiation has a more or 
less uniform distribution with a value of about .4 cal/cm^/min 
over North Africa, whereas the outgoing longwave radiation has 
^ strong gradient, starting with a minimum over the equatorward 
side of north Africa associated with high level clouds over the 
ITCZ and increasing northward with a maximum on the northernmost 
part of Africa associated with high surface tanperature in a 
clear area. Thus, a deficit is formed over the northernmost 
part of Africa. The outgoing longwave radiation in the model 
does not show a sharp northward gradient over this area and 
the values over this area are not large enough to produce a 
deficit. Among the possible reasons are the model-generated 
cloudiness associated with ITCZ and monsoon syste?i extends too 
far irarth and toat clouds are assumed to be black emitters and 
to fill an entire grid area as mentioned above. 

More clearly than for winter, some of the differences 
in the detailed structure of the pattern may be due to the 
failure of the model to generate low stratus and strato-cimulus 
clouds over the vest coasts of North Anerica, South Anerica, 
and Africa. Fbr example, instead of a local minimum, the model 
generates a local maximim over the west coast of North America. 
Ihe model-generated local maxima with values larger than 
.12 cal/cm 2 /inin extend too far south and extend into the 
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Southern Hemisphere over both the eastern Pacific and eastern 
Atlantic. Moreover, the model-generated cloudiness associated 
with the ITCZ system over the Atlantic extends too far 
These effects can be seen from Figs. 5b and 6b. Fig. 5b indi- 
cates cates that the local maxima with values larger than 
.45 cal/cm^/min of the absorbed solar radiation extend too far 
south into the eastern south Pacific. Fig. 6b indicates that 
the simulated local maxima outgoing longwave radiation over the 
south Atlantic is too far south in comparison with ^servations. 
The local mc«ima with values larger than .12 cal/can /min over 
the Northern Hemisphere also extend too far north (to 40N) over 
the eastern Pacific and eastern Atlantic, because the model- 
generated subtropical minimun cloudiness syston extends too far 
north. In other words, cyclonic activities are shifted too far 
north in the model. 

The model does not simulate a deep northwestward dip north 
of New Guinea, because the cloud systems generated by the model 
are shifted boo far east and have broader structures. The ob- 
served zero-balance isoline on the winter hemisphere is around 
lOS with northward dips over the east and west coasts of both 
South America and south Africa. The dips over the west cMst 
may be related to the stratus clouds and cloudiness associated 
with the ITCZ over these regions. The dips over the east 
coasts may also be related to the clouds over that area. The 
simulated zero-balance isoline is around 28S. The difference in 
solar radiation accounts for about an 8® difference in latitudes, 
but this still leaves a residue difference of about 10®. 

As far as the pattern is concerned, the maximum outgoing 
longwave radiation is confined between lOS and 30S, because the 
model-generated subtropical minimum cloudiness is located 
around lOS and 30S. There seams to be a coincidei^ of the zero- 
balance isoline with the poleward side of the max mum outgoing 
longwave radiation for both winter and sunmer, which ap^ars in 
both simulated and observed distributions. As for the inten- 
sity, we did not take the cloud fraction and cloud trananittance 
into account in our longwave radiation parameterization, with 
the result that we tended to underestimate the outgoing longwave 
radiation. 

On the whole, the differences in the detailed structure 
of the patterns and intensities are either due to; The model 
failinq to generate low-level stratus clouds over the west 
coasts of North America, South America, and Africa; or differ- 
ences in detailed structure of the cloud distribution ^ttems; 
or the assumption that clouds are black in infrared emission 
and cibsorption spectra, and fill an entire grid area. Of 
course, some of the differences are due to the natural varia 
bility of the atmosphere. The above comparison is based on an 
assirption that values as derived from Nimbus 3 measurement 
are correct. 
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Zonally-Averaqed Distributions and Global Mean 

TJie zonally-averaged distributions of net heat surplus 
and deficit with latitude, vAiich are important for the meridi- 
onal energy transport, are plottted in Fig. 7 for winter com- 
parison. Ihe values as derived from satellite measurement are 
for January; the values as simulated by model are for February 
It should be noted that the observations are taken near 
11:30 a.m. (9 a.m.) and 11:30 p.m. (9 p.m.) local standard 
time for Nimbus 3 (NOAA SR) measurements. Because the model 
systematically underestimates outgoing longwave relation, the 
net radiation balance is overestimated in oor^rison with re- 
sults of NCAA SR measurements, but it is overestimated only in 
the winter hemisphere in comparison with Nimbus 3 measuranents 



Fig. 7. Zonally-averaged radiation balance. 
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The global radiaticffi balance is tabulated in Table 1. 
Because the nxDdel-generated longwave radiation is systematically 
lower than values derived from Nimbus 3 measurenents, the model 
produces a large unbalanced radiation. Vfe have noted above 
that the overestimated cloud amount and cloud emissivity might 
account for much of the low outgoing longwave radiation. In 
addition^ the simulated tei^rature is colder than that of 
observation; for example, it is almost 10°C to 20®C too cold 
in polar regions vAiich will also tend to underestimate outgoing 
longwave radiation. 

Table 1. Global radiation budget of the earth-atmosphere system 


Outgoing 

Longwave radiation Radiation balance 

Albedo cai CB*^Bin“^ cal cb'^bIo"^ 


Observed 
(Jan. 21- 


February 3, 1970)% 

.501 

.361 

.283 

.337 

.024 

Slnulated 

{ February) 

.513 

.356 

.305 

.295 

.061 

obbcrved 

(July 16- 
July 31, 1969) 

.472 

.339 

.261 

.354 

-.015 

Simulated 

( Aunuiit ) 

.488 

.334 

.315 

.294 

.040 


SoUr radiation 
cal ca'^ln"*^ 
Incoalng Abaorbed 


Since the sea surface tanperature is fixed in the model, 
we suspect that most of the excess of the unbalanced radiation 
is lost to the oceans. This can be seen from Table 2, in which 
the heat balance of the ice- free ocean surface is tabulated. 

In the model, the sea surface temperature does not change 
whether the heat balance at the ocean surface is positive or 
negative. Ihe sensible heat flux, moisture flux, and net long- 
wave radiation flux in the model depend on the prescibed sea 
surface tanperature. Calculations, results of which are shown 
in Table 2, indicate that heat loss due to the above fl^es 
is less than heat gain due to solar flux by .083 cal/cm^/nin 
(.064 cal/on^/min) for February (August) simulations. Hie 
global mean, i.e., the above Vcilues weighted by fraction of the 
ice- free ocean (about 66%) is compatible with values of the un- 
balcinoed radiation at the top of the atmosphere. 
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Table 2. 


Heat balance of the ice-free ocean surface 
(cal cm"^min“^) . 



Net solar 
radiation 
absorbed 
by oceans 

Net loss 
through 
longwave 
radiation 

U>aa 

through 
sensible 
heat flux 

Loaa of heat 
through 
evaporat Ion 

Nat 

heating 
of ocean 

Global 
Ice- free 
ocean surface 
balance 

February 

.275 

.073 

.033 

.066 

.083 

.055 

August 

.168 

.042 

.018 

.044 

.064 

.042 


On an overall basis, the radiation balance of the earth- 
atmosphere systen simulated by using the GLAS GCM is ccnipatible 
with results as derived from Nimbus 3 observations in geograp^ 
ical distribution of the patterns but with differences in de- 
tailed structure of the patterns and in intensities. Vfe think 
we know seme of the causes of these discrepancies and are in 
the process of trying to eliminate them. 
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